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Abstract 

We report on preliminary measurements of time-dependent CP-violation asymmetries in charmless 
neutral B meson decays to K + K~K° (including resonant decays <pK° and / (980)K ), r]'K , ir°K®, 
KgK^Kg, Kg, p°K®, ujKg . The results are obtained from a data sample of up to 347 million BB 
pairs produced by e + e _ annihilation at the T(45) resonance collected with the BABAR detector at 
the PEP-II asymmetric-energy B-meson Factory at SLAC. 
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1 Introduction 



Measurements of time-dependent CP asymmetries in B° meson decays through a dominant Cabibbo- 
Kobayashi-Maskawa (CKM) favored b — > ccs amplitude have provided a crucial test of the mech- 
anism of CP violation in the Standard Model (SM) [2]. For such decays the interference between 
this amplitude and B°B° mixing is dominated by the single phase /3 = arg (—V c dV* b /VtdV t * b ) of the 
CKM mixing matrix. The BABAR measurement of (3 for these modes is sin2/3 = 0.710 ± 0.039 3 . 

In the SM, decays of B° mesons to charmless hadronic final states, such as 4K , /o(980)/f , 
K + K~K°, rj'K , tt°K®, KgKgKg, p°K®, ioKg, proceed mostly via a single loop (penguin) amplitude 
with the same weak phase as the b — > ccs transition 4 . In these modes, assuming the penguin 
dominance of b —* s transition and neglecting CKM-suppressed amplitudes, the time-dependent 
CP-violation parameter S (defined in Eq. ^ below) is expected to be sin2/3. However, CKM- 
suppressed amplitudes and the color-suppressed tree-level diagram, present in not pure-penguin 
modes (like f (980)K , K + K~K°, rj'K , ir°K®, p°K®, loK° s ), introduce additional weak phases whose 
contribution may not be negligible H3 IB] • As a consequence, only an effective S = sin2/3 e gf is 
determined. The deviation AS = S — sin2/3 has been estimated in several theoretical approaches: 
QCD factorization (QCDF) QCDF with modeled rescattering [TO, Soft Collinear Effective 

Theory (SCET) and SU(3) symmetry [El El El- The estimates are channel-dependent. QCDF 
and SCET models estimate AS to be positive in the most of modes. SU(3) symmetry provides 
unsigned bounds of the order \AS\ < 0.05 in the best case. 

Due to the large virtual mass scales occurring in the penguin loops, the possible presence 
of additional diagrams with new heavy particles in the loop and new CP-violating phases may 
contribute to the decay amplitudes. In this case the measurements of significantly larger AS are 
a sensitive probe for physics beyond the SM |Sj. Due to the different non-perturbative strong- 
interaction properties of the various penguin decays, the effect of New Physics (NP) is expected to 
be channel-dependent. 

In the SM, the decay B° — ► K®K® is expected to be dominated by the penguin b — ► d transition, 
and is potentially sensitive to the presence of NP in a way analogous to b — > s decays |13j . Neglecting 
CKM-suppressed amplitudes, the time-dependent CP- violating asymmetry parameters in this mode 
are expected to vanish (i.e. S = 0), while contributions from lighter quarks or supersymmetric 
particles could induce observable asymmetries [Hj . 

In this summary we report preliminary measurements of CP-violating asymmetries in all of the 
above mentioned loop-dominated decays. The data sample used consists of 347 million BB pairs 
(227 million for p K®,), recorded at the T(4S) resonance (center-of-mass energy \fs = 10.58 GeV). 
The data were collected with the BABAR detector 15 a at the PEP-II asymmetric-energy e + e~ 
collider. Detailed description for each analysis presented here are given in Refs. |16| and ^7]. In 
Ref. |15j we describe the silicon vertex tracker (SVT) and drift chamber used for track and vertex 
reconstruction, the Cherenkov detector (DIRC), the electromagnetic calorimeter (EMC), and the 
instrumented flux return (IFR). 

We reconstruct a B° decaying into the CP eigenstates K + K~K°, r/'K°, tt°K®, K g s K q s K g s , K° S K° S1 
p°Kg, ujK® (Bcp). From the remaining particles in the event we also reconstruct the decay vertex 
of the other B meson (B tag ) and identify its flavor. The difference At = tcp — t tSug of the proper 
decay times tcp and itag of the CP and tag B mesons, respectively, is obtained from the measured 
distance between the Bqp and B tag decay vertices and from the boost (/?7 = 0.56) of the e + e~ 
system. Due to the K® lifetime, the At for the modes tt°K®, KgKgK®, K®K® is obtained reliably 
by exploiting the knowledge of the average interaction point and from a global constrained fit to 



the entire T(45) — ► B°B° decay tree, including the constraint from the lifetime of the B° meson. 
The At distribution is given by: 

e -|At|/T 

F(At) = — [1 tAui±(1- 2w)(-ri f Ssm(Am d At) - C cos(Am d At))], (1) 

where rjf is the CP eigenvalue of the final state /, the upper (lower) sign denotes a decay accom- 
panied by a P° (P°) tag, r is the mean P° lifetime, Am d is the B°B° mixing frequency, and the 
mistag parameters w and Aw are the average and difference, respectively, of the probabilities that 
a true B° is incorrectly tagged as a P° or vice versa. The tagging algorithm has six mutually 
exclusive tagging categories ,3\ A non-zero value of the parameter C would indicate direct CP 
violation. 

2 Analysis Method 

Considering the K® candidates reconstructed in 7r + 7r~, we reconstruct the signal candidates com- 
bining Kg and K + K~, rj , it , K®K®, p°, lo, or another K® candidate. We also reconstruct the K® 
candidates in 7r°7r°, and they are combined with K + K~, rj , or K®K® (both in 7r + 7r~). Finally, 
we reconstruct the modes P° — ► K + K~K® and P° — » rj'K®, where a K® candidate is identified 
either as an unassociated cluster of energy in the EMC or as a cluster of hits in the IFR. The rj 
candidates are reconstructed in p°7 and in r/7r + 7r~, with r\ candidates in 7r + 7r~7r° (not considered 
in the signal candidates with Kg — > and K®) and in 77. 

We use the informations from the tracking system, the EMC and the DIRC to identify pions and 
kaons in the final state. Two kinematic variables are used to discriminate between signal decays 
and combinatorial background. The first is AE, the difference between the center-of-mass (CM) 
energy of the B candidate and the CM beam energy. The second is the beam-energy-substituted 

mass rriES = \J {s/2 + p ■ pb) 2 /Eq — p 2 B , where the B candidate momentum ps and the four- 
momentum of the initial T(4S) state (Eq,pq) are defined in the laboratory frame. In the 71"°-/^ 
and KgKgK® analyses, these variables are replaced by the invariant mass of the reconstructed B 
meson, tub, and the missing mass m m i ss = \q e + e - — qs\, where q e + e - is the four- momentum of the 
e + e~ system and qs is the four-momentum of the B candidate after applying a S°-mass constraint. 

Background events arise primarily from random combinations of particles in continuum e + e~ — > 
qq events (q = u,d,s,c). We reduce these with requirements on shape-event variables, like the 
angle 9t between the thrust axis of the B candidate in the T(4S) frame and that of the rest of 
the charged tracks and neutral calorimeter clusters in the event. In the fit we discriminate further 
against qq background with a Fisher discriminant J 7 01 a neural network (NN) which combines 
several variables that characterize the production dynamics and energy flow in the event ^H]- We 
study the background from other B decays using Monte Carlo (MC) simulated events. We take 
care of this background adding specific components in the fit. 

We obtain the CP-violation parameters and signal yields for each mode from extended maximum 
likelihood fits with the input observables AE, rriES; 3~ or NN, At (all modes) as well as the resonance 
mass and decay angle (p^K® and u)K®). In the case of modes with K\ the B mass kinematic 
constraint is necessary to determine the K® momentum so that tues cannot be exploited in the fit. 
In the fits, the likelihood for a given event is the sum of the signal, continuum and the P-background 
likelihoods, weighed by their respective event yields. 

In the K + K~K° analysis, we use an angular moment analysis to extract strengths of the partial 
waves in K + K~ mass bins. In this approach we rely only on the assumption that the two lowest 



Table 1: Preliminary fit results for each penguin mode. The second column gives /3 e g- or S, the 
third the result for the Aqp or C. The first errors given are statistical and the second systematic. 
See text for the description of <fiK° and fo(980)K° results. 



Mode /3 eff 


A C p 


K+K~K U 




0.361±0.079±0.037 


-0.034±0.079±0.025 






0.06±0.16±0.05 


-0.18±0.20±0.10 


/o(980)K° 




0.18±0.19±0.04 


0.45±0.28±0.10 



Mode 


s 


c 


V 'K U 








0.55±0.11±0.02 


-0.15±0.07±0.03 










0.33±0.26±0.04 


0.20±0.16±0.03 




K° s 






0.66±0.26±0.08 


-0.14±0.22±0.05 


P°K° S 








0.17±0.52±0.26 


0.64±0.41±0.25 








0.62lg;^±0.02 


— 0.43±g;||db0.03 










1 no+0.80 +0.11 
i.^o_ .73 _o.l6 


-0.40±0.41±0.06 



partial waves {S- and P-wave) are present, but make no other assumption on the decay model. 
Furthermore, for this mode we develop a novel technique based on a time-dependent Dalitz plot 
analysis, to take into account the variation of CP content and interference naturally in the fit. 
We use an isobar model where we include the resonances /o(980), 0(1020), Xq(1550), and XcO- In 
addition to resonant decays, we include non-resonant amplitudes. We extract /3 e g and Aqp = —C 
from the asymmetries in amplitudes and phases between B° and B° decays across the Dalitz plot. 

3 Results 

The fit results for the CP parameters are given in Tabled In K + K~K° analysis, we find that the 
trigonometric reflection at ir/2 — /3 c g is disfavored at 4.6cr, which is the first such measurement in 
penguin decays. From the angular moment analysis we find the P-wave fraction to be 0.29 db 0.03 
averaged over the Dalitz plot, and 0.89 ± 0.01 over the 4> resonance region (1.0045 < m K + K - < 
1.0345 GeV/c 2 ). We perform also a fit to \ow-K + K~ mass region (m K + K - < 1.1 GeV/c 2 ) in 
order to measure CP-asymmetry parameters for (j)K° and 70(980)^° components. This is the first 
measurement of the CP parameters in /o(980)-ftT° with /o(980) — ► K + K~ . 

All measurements reported here are statistics limited. The main sources of systematic errors 
are the Dalitz plot model (K + K~K°), the CP content of BB background, the reference shape 
modeling, the interference with other resonances (p°K®), and the SVT alignment (n°Kg, K®KgK®, 



The individual sin2/? e fj results are in agreement with the charmonium value. The sin2f3 e g for 
r]'K° is inconsistent with zero by 4.9 standard deviations. The S value for K®K® is consistent with 
zero. None of the modes studied exhibits non-zero direct CP violation. 



4 Conclusion 

We have presented preliminary results on mixing-induced and direct CP violation for loop-dominated 
B decays. We have studied several modes, reconstructing several sub-decays for each mode in order 
to increase the statistics. We have improved the analyses using better techniche, like time-dependent 
Dalitz plot analysis. 

All individual penguin modes are in agreement with the SM expectation. However, the sin2/? e fj 
measurements are lower than the charmonium value so that their average appears to be rather low. 
With the increase of statistics, we are very close to measure the CP violation in the penguin modes 
with b — > s transition. 
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